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ABSTRACT 

We investigate the jet launching process from accretion disks extending our recent study (paper 
I) to a truly bipolar setup. We perform axisymmetric MHD simulations of the disk-jet interaction 
on a computational domain covering both hemispheres, in particular addressing the question of an 
intrinsically asymmetric origin of jet / counter jet systems. Treating both hemispheres simultaneously, 
we overcome the equatorial plane symmetry boundary condition used in most previous studies which 
naturally fosters a symmetric evolution. For the magnetic diffusivity prescription we apply an a- 
parametrisation, considering both, globally models of diffusivity, and local models. We first approve the 
quality of our numerical setup by generating perfectly symmetric jets, lasting over a 1000s of dynamical 
time scales. We then disturb the hemispheric symmetry in the disk, and investigate the subsequent 
evolution of the outflow. The evolution first leads to a substantial disk warping with electric currents 
intersecting the equatorial plane. We investigate two models, i) a disk with (initially) different thermal 
scale height in both hemispheres, and ii) a symmetric disk into which a local disturbance is injected 
in one hemisphere. In both cases the disk asymmetry results in asymmetric outflows with mass fluxes 
differing by 10-20%. For case i) the jets are launched asymmetric right from the beginning. However, 
due to the symmetric diffusivity profile and restoring gravitational forces, the outflow asymmetry dies 
out after several hundred dynamical time scales. For case ii) the onset of asymmetry in the outflows is 
delayed as the disturbance injected at a larger radius first has to propagate to smaller radii where the 
bulk of the outflow is launched. For a local prescription of magnetic diffusivity (following mainly the 
local sound speed) the outflow asymmetry is persistent for long time, implying that diffusive effects 
play a major role when re-establishing hemispheric symmetric by gravitational forces. We find up to 
30% difference in mass flux between jet and counter jet for this setup, lasting over 1000s of dynamical 
time scales (i.e. lasting for the whole simulation). For comparison we have run simulations where jets 
are launched from a symmetric disk into an asymmetric coronal hemispheres of different density. As 
expected, the initial jet and counter jet are launched asymmetric, however, as soon as the jets have 
left the coronal asymmetry (within a few hundred dynamical time steps) a symmetric flow structure is 
established. In summary, our results suggest that the jet asymmetries in protostellar and extragalactic 
jets can indeed be generated intrinsically and maintained over long time by disk asymmetries and the 
standard jet launching mechanism. 

Subject headings: accretion, accretion disks - ISM: jets and outflows - magnetohydrodynamics - 
galaxies: active - galaxies: jets - stars: pre-main sequence 



1. INTRODUCTION 

Astrophysical jets as highly collimated beams of high 
velocity material and outflows of less degree of collima- 
tion and lower speed are an ubiquitous phenomenon. The 
jets are observed over a wide range of luminosity and spa- 
tial scale, originating from young stars, micro-quasars, 
or active galactic nuclei. The common understanding is 
that magnetohydrodynamics (MHD) are responsible for 
launching, accelerating, and collimating these jets. Fur- 
thermore, it is clear that accretion and ejection is related 
to each other - one efficient way to remove angular mo- 
mentum from a disk is to eject it vertically into a jet or 
outflow. 

Observational data have confirmed the co-existence of 
bipolar jets in most jet forming regions. Jet and counter 
jet appear, however, typically asymmetric in shape with 
very few exceptions (see e.g. Mundt et al. 1990 Ray 
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et al |1996[|Ray et al.|2007[ ). |Ellerbroek et al.|fl2013| ) find 
similar mean velocities in b oth lobes of the jet s observed , 
while they report asymmetric variations in mass outflow 
rates and velocities, and suggest that the jet launching 
mechanism on either side of the disk is not synchronized. 

One exception is the jet source HH 212 whic h ejects an 
almost per fectly symmetric bipolar structure ( jZinnecker 
et al.|1998 ), suggesting that the causal origin ol jet knots 



is located close to the central engine. On the other hand, 
if jets form naturally asymmetric, and if thus asymmetric 
jets would need special conditions to be formed, we may 
ask what is this kind of "natural" ejection process and 
what are the additional conditions for symmetry? 
A small number o f paper s have a addressed this topic. 



Chagelishvili et al. (1996) discuss the MHD origin of 
the one-sidedness as probably caused by a superposi- 
tion of a quadrupolar disk dynamo and a stellar dipole. 
The first bipolar jet simu lation has been published by 
von Rekowski et al. ( |2003| ) and follow-up papers. Re- 
cent simulations consider asymmetric ejections of stellar 
wind components from an offset multipole stellar mag- 
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netosphere ( Lovelace et al.pOlO ). Velocity as ymmetries 
in protostell ar jets have been investigated by Matsakos 
et al. ( |2012| ), in particular comparing intrinsic and ex- 
trmsic mechanisms, such as e.g. the field alignment be- 
tween stellar and disk field or the pressure distribution in 
the ambient cloud. Both effe cts seem to play a role and 
work on different time scales. Sauty et al. ( 2Q12[ ) discuss 
semi-analytical models indicating a counter-rotation of 
the inner jet and compare them to apply MHD simula- 
tions. 

To our knowledge, only very few numerical simulations 
investigating the bipolar launching of disk jets have been 
performed. It is therefore interesting to investigate the 
evolution of both hemispheres of a jet-disk system in or- 
der to see whether and how a global asymmetry in the 
large-scale outflow can be governe d by the disk evolu- 
tion. In a previous paper (paper I, |Sheikhnezami et al. 
2012 ) we have presented detailed model simulations lor 
the jet launching process, investigating a comparative set 
of various launching parameters such as magnetization 
(or plasma-beta), disk diffusivity and diffusivity scale 
height. In this paper (paper II) we continue our investi- 
gation concentrating on the bipolar character of jets and 
outflows, in particular on the fact that both jet compo- 
nents are mostly observed as asymmetric. We investigate 
the launching of bipolar outflows from initially symmet- 
ric conditions, and also from a disk structure with slight 
asymmetries. 

With launching, we denote the process which conveys 
material from radial accretion into vertical ejection, lift- 
ing it from the disk plane into the corona, thus estab- 
lishing a disk wind. Referring to our notation in paper I 
we like clarify again that with formation we denote the 
process of accelerating and collimating an already exist- 
ing slow disk wind or stellar wind in to a jet beam. This 
paper deals with both processes - launching and forma- 
tion. 

The present paper is organized as follows. Section 2 
describes the numerical setup, the initial and boundary 
conditions of our simulations with particular emphasis 
to the approach of simulating both hemispheres. The 
general evolution of jet launching is presented in section 
3, where we also present the simulation of a perfectly 
symmetric jet. Section 4 is then devoted to asymmetric 
jets and a parameter study comparing jets from different 
disk setups. We will discuss their symmetry properties 
and how symmetry can be broken by the intrinsic disk 
evolution. 

2. MODEL SETUP 

As discussed in detail in paper I, we model the launch- 
ing of MHD jets from slightly sub-Keplerian disks, which 
are initially in pressure equilibrium with a non-rotating 
corona. The main goals of paper I were i) to determine 
the mass ejection to accretion rate fraction for a vari- 
ety of disk physical characteristics such as plasma-beta, 
magnetic diffusivity, or disk scale height, and ii) similarly 
for the angular momentum flux, and iii) to determine the 
main geometry for these jets, such as the asymptotic jet 
radius and opening angle, or the size of the main jet 
launching area of the disk. 

In the present paper we extend our approach to simu- 
lations on a computational domain including both hemi- 
spheres. This allows us to investigate the truly bipolar 



launching and how the launching symmetry can affect 
the symmetry characteristics of jet and counter- jet. 
We apply the MHD code PLUTO (|Mignone et al.|2QQ7 



2Q12|, solving the time-dependent resistive MHD equa- 



tions as described in paper I. Again, the simulations are 
performed in a 2.5-dimensional setup in cylindrical coor- 
dinates (thus applying 3D axisymmetry ) . 



3. NUMERICAL SETUP - INITIAL , 
CONDITIONS 



BOUNDARY 



The major extension from paper I is that we now treat 
the truly bipolar launching of outflows. In general we 
apply the same initial disk structure and boundary con- 
ditions as before, however we extend the disk-outflow sys- 
tem across the equatorial plane into both hemispheres. 
Above and below the disk, respectively, a hydrostatic 
corona is prescribed in pressure balance with the disk 
gas pressure (thus, implying a density jump and entropy 
jump from disk to corona, see paper I). 

We apply a uniform grid with across the midplane and 
attached to that scaled grids, We run simulations ap- 
plying a different grid resolution (see Tab{T]). With the 
highest-resolution grid we resolve the disk scale height at 
the innermost radius with up to ~ 10 grid cells. 

3.1. Boundary conditions 

The main goal of this paper is to investigate the sym- 
metry of bipolar jets launched from a diffusive accretion 
disk. It is therefore essential to carefully check our nu- 
merical setup, in particular the internal boundary condi- 
tions describing the sink, in order to prevent numerical 
artifacts generating asymmetry. 

Compared to paper I, obviously the equatorial-plane 
boundary condition is now omitted. The disk itself may 
now evolve into a warped structure, breaking the intrin- 
sic hemispheric symmetry of disk and outflow. A disk 
mid-plane, if it exists, will not be necessarily along the 
equatorial plane. Further consequences are e.g. that the 
electric currents can now flow across the disk midplane. 

As discussed in paper I, for the outer disk radius we 
apply n outflow boundary condition. We feed the in- 
ner jet launching area by accretion from the inner disk 
areas. Since compared to paper I the physical extent of 
the computational domain is somewhat reduced, also the 
mass reservoir for disk accretion is smaller, which limits 
the disk mass evolution to comparatively shorter time 
scales. 

The remaining boundary conditions are equivalent the 
conditions similar we have applied in paper I, i.e. the 
outflow boundary conditions ( modified from the original 
code, see Porth fc Fendt|20 10) for the radial and vertical 
outer boundaries, and the axisymmetry boundary condi- 
tion along the rotation axis. The sink boundary condi- 
tions allow to absorb the mass and angular momentum 
of the accreting material. 

3.2. Initial conditions 

We apply the same basic initial conditions as in pa- 
per I. In addition, a few extensions to this setup were 
made in order to break the symmetry and to govern an 
asymmetric evolution. 

As in paper I, the initial magnetic field is prescribed 
by the magnetic flux function ip following |Zanni et al.| 
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Table 1 

Grid resolution of various parameter runs. The equidistant grid usually covers the highly diffusive regions of disk and disk wind, while the 
stretched grid covers the weakly diffusive, almost ideal MHD regions further distant from the launching region. Simulations discussed in 

this paper are labeled with a star '*'. 



run 


equidistant subgrid 


stretched su 


bgrids 






n r 


grid r-extension 


n r 


grids r-extension 




n z 


grid z-extension 


n z 


grids z-extension 


sbl *-sb4* 


800 


0.0 < r < 60.0 










1214 


-40.0 < z < 40 


oUU, oUU 


-100.0 < z < -40., 


40.0 < z < 100. 


sb5-sb7 


800 


0.0 < r < 60.0 










1214 


-40.0 < z < 40 


800, 800 


-100.0 < z < -40., 


40.0 < z < 100. 


cfl-cf5 


800 


0.0 < r < 60.0 










1214 


-40.0 < z < 40 


800, 800 


-100.0 < z < -40, 


40.0 < z < 100 


cfB 


1600 


0.0 < r < 60.0 










2428 


-40.0 < z < 40 


1200, 1200 


-100.0 < z < -40, 


40.0 < z < 100 


cf8 


2428 


0.0 < r < 40.0 










2428 


-10.0 < z < 10 


1000, 1000 


-50.0 < z < -10, 


10.0 < z < 50 


cfl4* 


1200 


0.0 < r < 50.0 










3600 


-60.0 < z < 60 








cfl6* 


2428 


0.0 < r < 40.0 










2428 


-20.0 < z < 20 


1000, 1000 


-50.0 < z < -10, 


10.0 < z < 50 















( 20071 ), 

q /r\ 3/4 rr? 5 / 4 

V lJ (m 2 + (z/r) 2 J 

where B z $ measures the vertical field strength at (r = 
H, z = 0). The (initial) field tension is determined by the 
parameter m. We apply m = 0.4 as in paper I. 

For the disk density and pressure we apply the same 
distribution as Eq. (6) of paper I. However, we also run 
models with an initially asymmetric disk structure, con- 
sidering a pressure scale height in the upper disk hemi- 
sphere e = e up = 0.15 different from the pressure scale 
height in the lower disk hemisphere e = edown = 

0.10. 

The initial disk corona follows the same distribution 
as in paper I. In order to compare intrinsic effects of 
asymmetric bipolar launching with external effects, we 
have, however, also run comparison simulations with an 
initially symmetric disk, but with a disk corona of differ- 
ent density /pressure in the upper and lower hemisphere, 
respectively. 

3.3. Units and normalization 

We apply the same code units and normalization as in 
paper I. Throughout this paper distances are expressed 
in units of the inner disk radius n, while and pd,i 
are the disk pressure and density at this radius, respec- 
tivel}^] For the sake of comparison to previous papers 
we may assume n = 0.1 AU for a protostellar jet and 
T[ = 10 Schwarzschild radii for an AGN jet. Velocities 
are measured in units of the Keplerian velocity ^k,o a t 
the inner disk radius. Thus, by assuming smaller inner 
disk radii, the outflow speed will become larger. Time is 
measured in units of t[ = ri/v^u which can be related 
to the Keplerian orbital period tk,[ = 27vt[. 

The (initial) disk aspect ratio e is the ratio of the 
isothermal sound speed to the Keplerian speed, both 
evaluated at disk mid plane, e = c s /vk- Pressure is given 

2 The index "i" refers to the value at the inner disk radius at 
the equatorial plane at time t = 



in units of pd,i = e 2 Pd,i^K i- The magnetic field is mea- 
sured in units of B x = B Z1 . We adopt i>K,i = 1, Pd,i = 1 
in code units. 

4. MAGNETIC DIFFUSIVITY 

In order to extend of the setup of paper I into a truly 
bipolar configuration we need to re-consider our model 
for the magnetic diffusivity rj(r,z). An obvious con- 
straint is that the prescription of diffusivity should not a 
priori influence the symmetry of the system. Since the 
disk structure may now evolve asymmetric in both hemi- 
spheres, the magnetic diffusivity must be able to follow 
such a disk evolution. Naturally, in asymmetric disks the 
disk midplane does not follow the equatorial plane. 

In general, we assume the diffusivity to be anomalous 
and of turbulent origin. Since we do not resolve the disk 
turbulence from first principles (e.g. by resolving the 
MRI), we apply a parameterized diffusivity distribution 
effectively following in principle an a-prescription. We 
neglect geometrically more complex distributions such 
as e.g. MRI-inactive dead zones of turbulent diffusivity 
or a more detailed local treatment of MRI turbulence 
jGanmiie 1996; Pessa h et al.|2007irG ressel 2010 ) [ Fromang 
|et al.||2U12[|Lesur et al.|2U13| |Bai & Stone||20rg| ~ 

As detailed in paper 1, we assume a diagonal diffusiv- 
ity tensor with the non-zero components = 7? p , and 
n rr = r] zz = 7/0, where we denote n p as the poloidal mag- 
netic diffusivity, and r]^ as the toroidal magnetic diffusiv- 
ity, respectively. The anisotropy parameter x — Vfi/Vp 
quantifies the different strength of diffusivity in poloidal 
and toroidal direction. Here, we apply x — 3.0 for all 
simulations. 

For the diffusivity profile, we have investigated several 
options which we will discuss in the following. Table [2] 
compares the parameter setups for the different magnetic 
diffusivity distributions applied in our simulations. 

4.1. A global magnetic diffusivity prescription 

A first option for the diffusivity distribution is to mir- 
ror the profile applied in paper I along the equatorial 
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Table 2 

Overview on our simulation runs, concerning the character of asymmetry imposed and the diffusivity distribution applied. Simulation 

runs discussed in this paper are labeled with a star '*'. 







Q\/"m "m f*i~v\r lnvf*$\ ki "n cr 

O V 111111*3 Ijl V kJL C-CXlVlll^ 


rj profile a 






sbl* 


reference run 


none (symmetric) 


h 






OLm,l 


= 3.0 


sb2* 


global disk asymmetry 


initial scale height, (e up ,edown) = (0.15,0.1) 


fi 










sb3* 


local disk asymmetry 


over pressure injected at r = 12 


Jl 






OLm,l 


= 3.0 


sb4* 


global disk asymmetry 


initial density contrast, (5 up , ^down) = (10 -3 , 10 -4 ) 


h 






&m,2 


= 3.0 


SDO 


local disk asymmetry 


over pressure injected at r = 18 


h 






OLm,2 


on 

— O.U 


sb6 


local disk asymmetry 


over pressure injected at r = 18 


h 






OLm,2 


= 1.0 


sb7 


local disk asymmetry 


over pressure injected at r = 18 


h 






%,2 


= 2.0 


cbl 


global disk asymmetry 


initial scale height, e up , edown — 0.15, 0.1 


h 










cb2-3 


global disk asymmetry 


initial scale height, e U p, edown — 0.15, 0.1 


/2, 


r 


1/3 


Ctm,2 


= 0.01 


cb4 


global disk asymmetry 


initial scale height, e up , edown — 0.15, 0.1 


/2, 


r 


1/3 


&m,2 


= 0.05 


cb5 


global disk asymmetry 


initial scale height, e up , edown — 0.15, 0.1 


h 






Ct m ,3 


= 0.05 


cb6 


global disk asymmetry, high res 


initial scale height, e up , edown — 0.15, 0.1 


/2, 


r 


1/3 


&m,2 


= 0.05 


cb8 


global disk asymmetry, highest res 


initial scale height, e up , edown — 0.15, 0.1 


/2, 


r = 


2/3 


&m,2 


= 0.01 


cbl4* 


global disk asymmetry 


initial scale height, e up , edown — 0.15, 0.1 


h 






Ct m ,3 


= 0.1 


cbl6* 


global disk asymmetry 


initial scale height, e up , edown — 0.15, 0.1 


h 






Ct m ,3 


= 0.1 



t=0 t=100 t=1000 t=3000 




Figure 1. Time evolution of the bipolar jet-disk structure for reference simulation sbl applying a fixed-in-time and fixed-in-space magnetic 
diffusivity distribution Eq. [2[ Shown is the evolution of mass density (color) and the poloidal magnetic field (contours of poloidal magnetic 
flux V(r, z)) for the flux levels V = 0.01, 0.03, 0.06, 0.1, 0.15, 0.2, 0.26, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95, 1.1, 1.3, 1.5, 1.7, for the dynamical 
times steps t = 0, 100, 1000, 2000, 3000. 

plane, 



r] p (r,z) = a m ,i/i(r,2:) = a mA v A ,oH exp ^-^J, (2) 

with the Alfven speed ^a,o = va(t,z = 0) and the 
disk thermal scale height Hq = H(r) = er = cs(r, z = 
0)/vk(^j z = 0) 7 measured at the midplane and at time 
zero. Essentially, the diffusivity profile Eq. [2] is geomet- 
rically tied to the equatorial plane, potentially inducing 
hemispheric symmetry of the disk and the outflow. 



As in paper I, we allow for a diffusivity scale height 
H v larger than the ther mal scale heigh t Hq (see discus- 



sion in paper I). In fact, Gressel| ([2010) who investigated 
the MRI-induced turbulence ol accretions disks by high- 
resolution box simulations, finds an increasing level of 
turbulence with disk height. His simulations attest a 
maximum level of turbulence at about 2-3 disk pressure 
scale heights, which are in nice agreemen t with our model 
appro a ch. More recent si mulations by Beck with et al. 
(2011); Simon et al. (2013) indicate similar scale heights 
for the turbulent stresses. 
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Figure 2. Time evolution of the mass fluxes for reference simula- 
tion sbl. Shown is evolution of accretion rate and the mass ejection 
rates from the upper (solid) and lower (dashed) disk surfaces (all 
in code units). Ejection rates are measured in the control volumes 
with n = 1.5 and r<i = 10.0, while the accretion rate is vertically 
integrated at r = 10. 

In Eq. [2j both va and cs can be chosen time- 
independent (as in p aper I), or evolving in time (as e.g. in 
Murphy et at1|2010| ). For the sake of comparison, in the 



present bipolar reference simulation (denoted by sbl) we 
apply a magnetic diffusivity profile /i(r, z) fixed in time. 

4.2. A local prescription of magnetic diffusivity 

In the present paper we investigate an asymmetric 
evolution of the disk-jet structure in both hemispheres. 
Therefore, we need to apply a new, truly bipolar setup 
in which the leading parameters for the diffusivity pro- 
file are not defined with respect to the equatorial plane 
anymore. It is clear that for an asymmetric disk evolu- 
tion, a hypothetical disk mid-plane will not coincide with 
the equatorial plane. Therefore, in order to consider the 
spatial evolution of the asymmetric disk structure, it is 
essential to apply a local prescription of diffusivity. 

Among other options, one possibility is to relate the 
magnetic diffusivity to the local pressure or density. One 
may assume that the turbulent Alfvenic pressure is pro- 
portional to the therm al pressure as discussed e.g. by 
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Ouyed fc Pudritz ( |l997 ). In particular, Fendt & Cemeljic 



UU2| ) suggested that if magnetic diffusivity is interre- 
lated to the turbulent Alfvenic pressure, the magnetic 
diffusivity profile will follow a power law r] m ~ p 1 / 3 . In 
this casse the diffusivity is proportional to the local sound 



2/2(^2) «m,2p r 



(3) 



The profile with r ~ 1/3 results in a relatively broad ver- 
tical diffusivity profile, wider than the previously used 
exponential profile Eq[2] It has been shown that this 
will impact the jet acceleration and collimation (Fendt & 
|Cemeljic|2QQ2[ ) , maybe more than it effects the launching 
mechanism itself. For comparison we have also applied 
different power laws, such as T = 2/3, or even steeper 
profiles. The drawback of the simple profile [3] is that 
the diffusivity decreases along the disk is decreasing with 
radius. The outer disk material is therefore strongly cou- 
pled to the magnetic field, leading to an super efficient 
angular momentum removal, rapid accretion, and, thus, 
a short life-time of the outer disk. On the other hand, 
in case of MRI driven turbulence, the MRI activity is 
expected to cease for large radii, where "dead zones" for 
active accretion may exis t ( Gammie|1996||Beckwith et al 
20TT| [SilnTm et al (2013b . 



"We have favored another option for the magnetic diffu- 
sivity which allows for a disk diffusivity increasing with 
radius, 



7? p (r, z) a m , 3 / 3 (r, z) a m , 3 # L (r, z) 



1 



1 1 Hjjj^z) 
^ H ' 



(4) 



Thus, we apply a density- weighted "local disk scale 
height" (r, z) mainly following the local sound speed 
in the gas flow, 



H L (r,z) = p(r,z) 




P(r, z) 



1 



p(r,z) v Kep (r) 



(5) 



and a quenching term in order to avoid numerically prob- 
lematic diffusivities. We typically choose Ho = 0.5, 
o-p = 3/2, cr r 3/2,or H 0.1, a p 3/2, a r 5/2 
(simulation run cbl4). 

Essentially, both prescriptions for the magnetic diffu- 
sivity Eqs. [3] and [4] allow for a smooth transition from 
accretion to ejection, and, also to follow the changes 
in the local disk structure. As the outflow density de- 
creases along the streamlines, also the outflow diffusiv- 
ity decreases. For low densities towards the asymptotic 
outflow, the ideal MHD will be approached. A physical 
motivation can be the following. Turbulently diffusive 
disk material is lifted from the disk into the disk, and is 
further accelerated along the outflow while, however, the 
turbulent motions decay. In paper I we have estimated 
the scale height where the turbulence will be damped to 
be several d isk thermal sc ale heights (see our discussion 
above citing |Gressel|2010| ) . 

5. TEST CASE - LAUNCHING SYMMETRIC BIPOLAR JETS 

We first present jet launching simulations resulting in 
symmetric bipolar jets. The first example is the evolu- 
tion of jets following a magnetic diffusivity description 
fixed in time and space (case sbl). This case serves as 
reference simulation for this paper, and also allows for a 
comparison to the one-hemispheric simulations of paper 
I. This inflow-outflow evolution is shown in Fig. [I] 

As we see, the outflow evolves perfectly symmetric in 
both hemispheres for almost 3000 dynamical time steps, 
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Figure 3. Time evolution if the bipolar jet-disk structure for simulation sb2 applying a fixed-in-time diffusivity profile Eq.[2l and evolving 
from an asymmetric initial state with different thermal scale heights for the upper and lower disk hemisphere, e up = 0.15 and Cdown — 0.1. 
We show the evolution for time t = 0, 100, 1000, 2000 of the mass density (colors) and the poloidal magnetic field (lines), i.e. contours of 
poloidal magnetic flux * with flux contours = 0.01, 0.03, 0.06, 0.1, 0.15, 0.2, 0.26, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95, 1.1, 1.3, 1.5, 1.7. 
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Figure 4. Time evolution of the mass fluxes for simulation sb2. Shown is evolution of accretion rate and the mass ejection rates from 
the upper (solid) and lower (dashed) disk surfaces (all in code units). Ejection rates (middle) are measured in the control volumes with 
7*1 = 1.5 and T2 = 10.0, while the accretion rate (top) is integrated at r = 10. For comparison we show the asymptotic (vertical) mass 
fluxes through the jet, integrated from r\ = 2.0 to r2 = 50.0 at z = —75, 75. 



until the outer disk starts to deviate from symmetry due 
to numerical effects. However, even for these late time 
steps, the inner disk, which is the main jet launching 
area, is still highly symmetric. 

In particular this could be seen in the time-evolution 
of the mass fluxes. Figure [2] shows the accretion rate 
and the ejection rates into the two hemispheres. The 
ejection rates are integrated from r = 2 to r = 20 at 
z = ±3iJ(r), while the accretion rate is integrated from 
z = — 3H till z = 3H at r = 10. It is not surprising 
that due to the symmetric diffusivity profile which re- 
mains fixed in time together with symmetric initial and 
boundary conditions, we find a symmetric evolution of 
both the disk and the outflows. The symmetric bipolar 
jet structure we obtain does mainly serve as a test case 
for the numerical setup, in particular the sink boundary 



conditions. While the symmetric evolution is expected 
on physical grounds, even a little numerical failure in 
the setup would lead to asymmetry quickly. 

In comparison to our one-hemispheric simulation in pa- 
per I (reference run easel), the mass fluxes we measure 
now are quite similar. As accretion rate for the bipolar 
simulation we measure M acc (t = 2000) ~ 0.03, while for 
the previous reference easel it was M acc (t = 2000) ~ 
0.015 The outflow rates for the bipolar simulation are 
M ejec (t = 2000) 



with the M e j ec (t 



t 0.01 in each direction, which agrees 
2000) = 0.007 in paper I nicely. How- 
ever, due to the smaller grid extension, the disk mass 
reservoir is smaller, leading to a faster decay of the disk 
mass (see the outer disk structure at t = 3000). There- 
fore, disk accretion and mass ejection decay faster as well. 
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Figure 5. Disk evolution for the jet-disk system sb2. Shown 
is the radial velocity v r (r, z) for the dynamical time steps t = 
40, 100, 1000, 2000 (from top) indicating the asymmetric evolution 
of the disk and the disk wind. 

6. BIPOLAR JETS FROM ASYMMETRIC DISKS 

Here we present simulations in which we have dis- 
turbed the internal hemispheric disk symmetry. We have 
applied two options for disturbing the disk symmetry - 
either by prescribing a global asymmetric initial state, or 
by injecting a localized overpressure (an explosion at cer- 
tain time). For both alternatives we apply a symmetric 
magnetic diffusivity prescription Eq. [2j constant in time. 

6.1. An asymmetric disk scale height 

Option I is to prescribe an initial disk structure with a 
global pressure asymmetry in both disk hemispheres. In 
our model, we achieve this by applying a different ther- 
mal disk scale height for the initial disk in each hemi- 
sphere. In simulation sb2 we have applied e = H/r = 0.1 
for the upper hemisphere and e = H/r = 0.15 for the 
lower hemisphere. Consequently, we have a density and a 
pressure jump across the equatorial plane, AP/P = 0.2. 



The disk turns into an asymmetric evolution right from 
the beginning, evolving into a warped structural along 
the midplane (Fig. |3|. A series of warps is visible along 
the disk with warp amplitudes of a few local disk scale 
heights. After about 1000 dynamical time steps the warp 
amplitudes start to decrease - first along the inner disk, 
while the outer disk still in a warped state. 

The disk asymmetry is reflected in the jet evolution. 
Along with the initial asymmetric disk evolution, the jets 
launched from the inner disk are asymmetric as well. 
This is clearly visible in the poloidal magnetic field struc- 
ture (Fig. [3]), but also in the mass fluxes we measure. 
We also note a different time scale in the jet propaga- 
tion. The upper jet reaches the grid boundary earlier 
than the lower jet which is delayed by about At = 10% 
(note that this obviously depends on the grid size, and 
happens rather early at t < 50). 

Figure [4] shows the time evolution of the mass fluxes. 
Comparing both jet fluxes we find that at early stages, 
t < 700, the lower outflow carries about 80% of the mass 
load of the upper outflow. From t ~ 700 — 1500 an 
asymmetric inflow- out flow system is established. The 
differences in mass flux are now about 10%. After t ~ 
1500, the variation in the outflow rates decreases, and 
the inner disk has established a symmetric structure. 

The same behavior is also visible in the velocity distri- 
bution. Figure [5] shows the complex velocity field of the 
inner disk. Accretion starts first in the innermost disk 
regions. The disk asymmetry is well reflected in the v r - 
distribution. However, after t = 1000 the system turns 
to a symmetric geometry. We also see the highest accre- 
tion velocities in the upper disk layers. One may indicate 
this as layered accretion, however, due to the higher den- 
sity in the lower disk layers, the radial mass flux is more 
almost distributed over the disk height. 

We believe that the disk evolution into a symmetric 
state is given by two reasons. Firstly by the restoring 
force of the symmetric gravitational potential, and sec- 
ondly by the symmetry and the time-independent pre- 
scription of the magnetic diffusivity. The latter is in- 
teresting since it is a second-order effect only, as the 
magnetic diffusivity does not provide a force term in 
the MHD equations which could directly re-configure the 
disk structure. The diffusion time scale r v = AZ 2 /n ~ 
H 2 /rj ~ is about 100 dynamical time scales for rj ~ 
0.01, H ~ 1 at radius r = 10. 

So far we have concentrated on the inflow-outflow dy- 
namics close to the jet launching region. We now consider 
the evolution further downstream the outflows. This is 
interesting as it is this part of the outflow which is in prin- 
ciple accessible by the observations. Figure JI] (bottom) 
shows the mass fluxes of jet and counter-jet far from the 
source, integrated from r = 2 to r = 50 at z = 75, and 
z = —75, respectively. We find that the mass flux asym- 
metry of about 5% at the launching region propagates 
to the asymptotic region, where we find a similar mass 
flux difference. The time lag between launching region 
and asymptotic domain of about At = 300 = 1200 — 900 
can be explained by the propagation period of the accel- 
erating outflow of velocities v z ~ 0.2 — 0.8 propagating a 
distance Az = 75. The maximum jet velocity is achieved 
along a narrow cone between r = 10 and r = 20 (at 

3 Note, however, the axisymmetric setup 
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Figure 6. Time evolution of the inner disk accretion for simula- 
tion sb3, applying the diffusivity distribution fixed in time Eq. [2] 
An overpressure is added at t = 400 in the upper disk hemi- 
sphere at radius r = 12, lasting for few rotations. Shown is 
the evolution of the mass density (color) for the dynamical times 
t = 400, 404, 420, 600, 2800. 



z = 100). The outflow velocity of the lower jet is about 
1.2, while for the upper jet it is 1.1 times the Keplerian 
speed at the inner disk radius. The bulk of the mass flux 
of the outflow is, however, located between r = 15 and 
r = 30, and is a factor 4 higher than for the fast flow and 
moves with a speed v z ~ 0.4. 



6.2. A localized disk asymmetry in one disk hemisphere 

We start simulation sb3 with a symmetric initial disk 
structure. As for reference case sbl the outflows evolve 
into a symmetric jet-counter jet structure (Fig. [7]). How- 
ever, at t = 400, when a quasi steady state of the inner 
inflow-outflow structure is reached, we disturb the sym- 
metry of the disk structure by inserting a localized over 
pressure in the upper disk hemisphere, probably similar 
to an explosion (see Fig. [6J. This injection is localized 
within a box of size (Ar x Az) = (1.5 x 0.4) located 
at (r, z) = (11.25,1.2), and is switched on for At = 20, 
corresponding to 0.1 of an orbital period at this radius. 
The injected material has on average a 20 times higher 
density and a 2000 times higher pressure compared to 
the ambient disk material. 

The injected material disturbs the disk symmetry as 
it expands across the disk. The disturbance is slowly 
advected into the jet launching region. However, we ob- 
serve that the disk asymmetry decays faster than it is 
advected along the disk into the jet launching area at 
small disk radii. This is easy to understand. The expan- 
sion is roughly with sound speed, while the advection 
happens a sub-sonic speed v r ~ e^Kep- Figure [6] shows 
the disk accretion velocity evolution. The expansion fol- 
lowing the "explosion" first penetrates the whole disk 
in vertical direction, before the over-density is advected 
slowly inwards. Thus, the symmetry of the inner disk is 
not really affected by the explosion in the upper hemi- 
sphere at r = 12. The dominant outflow from the inner 
disk is only weakly affected. When the size of the launch- 
ing area reaches the explosion site, the asymmetry has 
almost disappeared. 

Up to t = 1000 the mass fluxes show an asymmetry 
on a 10%-level. Also a slight structural asymmetry is 
visible (see the 6 th flux surface contour in Fig. [7| which 
has propagated from the disk surface to a distance z ~ 70 
at this time. 

That part of the outflow which originates from the ini- 
tially asymmetric part of the disk (where the blob was 
injected), starts indeed asymmetric. As the mass flux 
launched from larger disk radii is small compared to the 
jet launched from the inner part, the total jet mass fluxes 
into the two hemispheres differ only marginally. Figure 
[8] show the time evolution of the mass fluxes. The ejec- 
tion rates are integrated along the disk surface between 
r = 1.5 and r = 10 thus inside the radius where the ex- 
plosion happens. Still the explosion disturbs the launch- 
ing site such that we see a 5-10% effect in the outflow 
mass fluxes. This effect is somewhat delayed from the 
explosion time since the disturbance needs time to be 
advected inwards. The asymptotic outflow rates, inte- 
grated from r = 2 to r = 50 at z = ±75, show only 
marginal differences, however. Note that the flux surface 
passing through (r, z) = (50, 75) anchors at r = 5 in the 
disk, and thus in a weekly disturbed region. 

In summary, although the initially symmetric disk is 
clearly disturbed by the asymmetric explosion, the asym- 
metry decays faster than it propagates to the inner jet 
launching site. Thus, the asymptotic, collimated jet, 
which originates in the inner disk area is only marginally 
asymmetric, moreover as the main mass flux is launched 
along the innermost field lines. 
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Figure 7. Time evolution of the bipolar jet-disk structure for a simulation sb3 with the fixed-in-time diffusivity distribution in time, 
Eq. [2] but with a localized overpressure added at t = 400 for At = 20. See Fig. [6] for a higher resolution of the inner disk area. We show 
the evolution for time t = 0,100,1000,2000,3000 of the mass density (color) ana the poloidal magnetic field (lines), i.e. contours of the 
magnetic flux, with flux levels ^ = 0.01, 0.03, 0.06, 0.1, 0.15, 0.2, 0.26, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, 0.95, 1.1, 1.3, 1.5, 1.7. 
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Figure 8. Time evolution of the mass fluxes for simulation sb3. Shown is evolution of accretion rate and the mass ejection rates from 
the upper (solid) and lower (dashed) disk surfaces (all in code units). Ejection rates (middle) are measured in the control volumes with 
7*1 = 1.5 and r2 = 10.0, while the accretion rate (top) is integrated at r = 5. For comparison we show the asymptotic (vertical) mass fluxes 
through the jet, integrated from r\ = 2.0 to r2 = 50.0 at z = —75, 75. 
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Figure 9. Time evolution of the inner disk structure for simulation run cbl4. Here the diffusivity profile is prescribed as a local Eq. |4] 
Shown is the evolution of magnetic diffusivity (color) and the poloidal magnetic field (contours of poloidal magnetic flux ty(r,z)) for the 
dynamical times steps t = 50, 500, 1000 (from left). 
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Figure 10. Time evolution of the jet mass flux rates for simulation 
cbl4. Shown is also the launching mass flux measured at 3 disk 
scale heights (top) and the asymptotic fluxes (bottom) integrated 
at z = ±50 from r = 2 to r = 40. 




7. 



A LOCAL MAGNETIC DIFFUSIVITY MODEL AND 
BIPOLAR JET LAUNCHING 



In order to allow the disk and jet to follow a truly 
asymmetric evolution, without the restoring effects of a 
symmetric, global magnetic diffusivity distribution, we 
have applied a local description of diffusivity which fol- 
lows the local evolution of disk and outflow (see Tab. [2]). 

The results of this section are somewhat preliminary, 
as a physically self-consistent parametrisation for a local 
magnetic (turbulent) diffusivity is not available. In or- 
der to investigate the main features of the local approach, 
we have applied two diffusivity distributions. One follows 
Eq. [3j the other one refers to Eq. [4j The main difficulty 
one faces with a local description is a feedback mecha- 
nism such that low densities lead to a lower diffusivity, 
thus stronger matter-field coupling, thus more efficient 
angular momentum removal and a faster accretion, which 
leads to even lower densities. 

Moreover, the local prescription Eq. [3] gives a diffusiv- 
ity profile decreasing with radius, contrary to the usual 
choice in the literature, Eq. 151 which is also realized in 
Eq.gj 

We find that the feedback mechanism is most efficient 
in simulations with a strong coupling between diffusivity 
and density, such as a diffusivity profile Eq. [3] In the ex- 
treme cases we see that the feedback mechanism may lead 
to a temporary gap opening over the innermost disk radii. 
As a result, the accretion in this area is strongly episodic, 
and, subsequently also the mass ejection. We will present 



a detailed investigation of this feedback mechanism in a 
forthcoming paper. 

Applying a more sophisticated prescription of the mag- 
netic diffusivity following Eq. HJ the over-all accretion- 
ejection evolution is more similar to the picture estab- 
lished by the simulations using a global diffusivity pro- 
file. The fundamental difference is, however, the longer 
lasting and more persistent asymmetry in the disk and 
the outflows. 

Figure [9] shows the time evolution of magnetic diffusiv- 
ity and the magnetic field for a simulation with a more 
sophisticated local setup of magnetic diffusivity following 
Eq. [4j We clearly see how the magnetic diffusivity fol- 
lows the structure of the disk and the outflow. The disk 
"warping" is seen diffusivity distribution as well. Mass 
loading and matter-field coupling depends on the local 
disk properties (defined by Eq. Hi). Since the diffusiv- 
ity profile is broader in vertical direction, mass loading 
and angular momentum removal is more efficient, estab- 
lishing higher accretion and outflow rates - in agreement 
with our results in paper I. 

The outflow mass fluxes develop a clearly asymmetric 
structure. Figure [10] show the mass fluxes integrated 
from r = 1.5tor = 10, along a surface z = 0.3r parallel 
to the initial disk surface. Jet and counter jet injection 
mass flow rates differ substantially - by about 30% over 
at least 1000 time steps. The asymptotic jet mass fluxes 
integrated from r = 2 to r = 40, along a surface z = ±50 
are somewhat lower (due to the fact that some of the 
material injected into the outflow leaves the grid in radial 
direction), and still differ on a 15%-level. 

8. IMPACT OF THE ENVIRONMENT 

This section is devoted to the question whether jet 
asymmetries are due to external of internal interaction. 
It seems clear that the conditions in the ambient gas the 
jet is penetrating, will certainly affect the expansion and 
the collimation on the large scales. 

Here we consider a model setup, in which the bipolar 
outflow is launched intrinsically symmetric from a (sym- 
metric) disk into an asymmetric star-disk corona. We 
do this by prescribing a different initial coronal density 
(and thus pressure), simply choosing S up = 10 -4 and 
5 up = 10~ 3 . 

Figure |TT| shows the time evolution of density and mag- 
netic fieldfor such a setup. 

As expected the initial outflow is highly asymmetric, 
more than the outflows disturbed by intrinsic disk asym- 
metries. However, the asymmetry is clearly transient. As 
soon as the initial coronal material is swept out of the 
computational domain, the accretion-ejection system re- 
turns to hemispheric symmetry. This is well visible in 
the mass flux evolution (Fig. 12), which shows a dras- 
tic asymmetry during the firstnundreds of time steps. 
Again it is interesting to see various time lags in the flow 
evolution. Firstly, the time lag between the launching 
time of the initial asymmetry and the arrival time of 
asymmetric features in the asymptotic region. Secondly, 
the time lag between the time step when the disk-outflow 
system has returned to symmetry, and the time when the 
outflow symmetry has reached the asymptotic jet. 

Note that in these simulations the disk structure itself 
stays rather symmetric, in difference to the simulations 
discussed above, where we see the disk evolving into a 
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warped structure. However, we can see back reaction 
of the asymmetric outflow onto the disk structure in this 
approach. The two asymmetric jets drive a different elec- 
tric current system in both hemispheres. Both the cur- 
rent systems are connected within the disk - subsequent 
MHD forces therefore slightly distort the symmetric hy- 
drodynamic disk structure. The deviation from symme- 
try is, however, not as strong as for the previous cases in 
which we start with an initially asymmetry disk. 

9. RESOLUTION STUDY 

Finally we present example results of our resolution 
study in brief. Figure [13] shows the inner part of the 
disk at the dynamical time step t = 130 when the inner 
disk evolution is rather violent. We compare simulation 
cbl4 applying our standard resolution of Ar = 0.0417, 
Az = 0.0333 in the disk, with simulation cbl6, applying 
a 3 times higher resolution, Ar = 0.0165, Az = 0.00824. 
Although the higher resolution run shows somewhat 
more sub-structure, such as internal shocks or a more 
peaked diffusivity distribution for r < 3, the main fea- 
tures of the disk dynamics are just the same. In partic- 
ular the disk height is similar, as well as the structure 
and opening angle of the magnetic flux surfaces. Also 
the radially structured features of the outflow close to 
the disk is very similar in both hemispheres. The over- 
all mass fluxes measured in the simulations are similar as 
well, however, the jet-counter jet asymmetry is somewhat 
more pronounced in the high-resolution run. 

We note that high-resolution simulations are particu- 
larly difficult to manage in diffusive MHD since the CF 
time stepping condition At v < (Ax) 2 /rj. 

10. CONCLUSIONS 

We have presented results of MHD simulations investi- 
gating the launching of jets and outflows from a magnet- 
ically diffusive disk in Keplerian rotation. The time evo- 
lution of the accretion disk structure is self-consistently 
taken into account. The simulations are performed in 
axisymmetry applying the MHD code PLUTO. 

Based, on paper I which studied how magnetic diffusiv- 
ity and magnetization affect the disk and outflow proper- 
ties, such as the mass and angular momentum fluxes, jet 
collimation, or jet radius, the present paper investigates 
the hemispheric symmetry of outflows. We have setup 
a numerical scheme that can treat the outflows from an 
equatorial disk in bipolar direction. The setup has been 
carefully checked against numerical artifacts triggering 
asymmetry in the disk-outflow evolution. 

We disturb the disk-outflow bipolar asymmetry apply- 
ing several approaches. Bipolar symmetry is re- achieved 
on the long term by gravitational forces, in particular 
the component vertical to the disk, and (in some cases) 
a symmetric magnetic diffusivity prescription. 

In particular, we have obtained the following results. 

(1) A test case with a symmetric setup gave a perfectly 
symmetric bipolar evolution of disk and outflow for sev- 
eral thousand rotations, until numerical effects by the 
outer disk boundary condition start to disturb the sym- 
metry of the outermost disk. The measured mass fluxes 
compare well with the one-hemispheric simulations in pa- 
per I. 

(2) We then applied various options to disturb the 
outflow bipolar symmetry. First we prescribed an ini- 



tially asymmetric disk applying a different pressure scale 
height in both disk hemispheres, e up = 0.1, ed OW n = 0.15. 
In general the disk structure evolves into a warped config- 
uration, with warp amplitudes of a few initial disk scale 
heights. Electric currents are driven across the equatorial 
plane. The mass fluxes of the resulting outflows differ by 
about 10% over less than 1000 time steps, until a sym- 
metric launching is achieved again. The outflow rates at 
larger distance from the source differ similarly, however, 
asymmetry last longer, just because of the propagation 
time of the material launched from the disk. 

(3) We then started the simulation with symmetric ini- 
tial setup, but disturbing the symmetric disk structure 
by a localized (in time and space) explosion at the time 
when a symmetric outflow is already well established. As 
for (2), a asymmetric outflow arises after the injection, 
however, after a time delay due to the propagation time 
of the asymmetric injection along the disk towards the 
main jet launching area. On the long term a symmetric 
outflow is re-established. 

(4) We then investigated the disk-jet evolution ap- 
plying a local prescription of magnetic diffusivity n = 
77 (r, z, t) following the local disk evolution. We find that 

- as discussed in paper I - the distribution of magnetic 
diffusivity affects the disk evolution (and subsequently 
the outflow evolution) essentially, as it governs the cou- 
pling between matter and field, and, thus the angular 
momentum evolution, the disk accretion, the magneto- 
centrifugal acceleration, and also the mass loading of the 
outflow. We first applied a simple power law diffusiv- 
ity distribution as motivated by earlier papers. We ob- 
served a strong feedback such that for low densities an 
almost ideal MHD situation leads to super-efficient ac- 
cretion with weak outflows and strongly decaying disk 
mass. These features deserve further attention which 
will be the focus of a forthcoming paper. 

(5) For a more sophisticated magnetic diffusivity distri- 
bution following the density-weighted local sound speed, 
and an increasing diffusivity with radius, the accretion 
evolution is more persistent. We were able to follow the 
accretion-ejection process over many thousands of dy- 
namical time steps. The most interesting result is that 
the bipolar asymmetry of jet and counter jet is long- 
lasting. We find that in this case the warped structure 
of the disk survive many dynamical time scales also in 
the inner disk. We interpret this as due to the lack of 
a symmetric diffusivity distribution. There are the same 
restoring forces of gravity as in the case of a symmetric 
diffusivity profile. However, the matter is more directly 
coupled to the distorted magnetic field which in principle 
opposes the return to symmetry. In the end we observe 
a persistent difference in the jet-counter jet mass fluxes 

- up to 30% and lasting longer than the simulation run 
time. 

(6) In order to compare between internal and external 
effects we also investigated the launch of a jet outflow 
from a symmetric disk into an asymmetric ambient gas 
distribution. As probably expected the initial outflow is 
strongly asymmetric, with 20% different mass fluxes for 
jet and counter jet, however as soon as the outflow has 
penetrated the asymmetric corona (i.e. when the outflow 
has left the grid, and the initial condition has swept out 
of the domain, the outflow returns to symmetry (after 
about 2000 dynamical time steps). In comparison to the 
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Figure 11. Time evolution of the bipolar jet-disk structure for simulation sb4. Here the outflow is launched in hemispheres with different 
density contrast. Shown is the evolution of mass density (color) and the poloidal magnetic field (contours of poloidal magnetic flux \P (r, z)) 
for the dynamical times steps t = 0, 100, 500, 1000, 3000. 
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Figure 12. Time evolution of the mass fluxes for simulation sb4. Shown is evolution of accretion rate and the mass ejection rates from 
the upper (solid) and lower (dashed) disk surfaces (all in code units). Ejection rates (middle) are measured in the control volumes with 
7*1 = 1.5 and r2 = 10.0, while the accretion rate (top) is integrated at r = 5. For comparison we show the asymptotic (vertical) mass fluxes 
through the jet, integrated from r\ = 1.5 to r2 = 40.0 at z = —80, 80. 



simulations discussed about, the initial disk structure is 
symmetric and does not evolve into a warped structure. 
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Figure 13. Resolution study for our setup. Subset of the in- 
nermost disk-jet area which is most critical to resolution issues. 
Shown is the magnetic diffusivity distribution following Eq. [l] for 
simulation run cbl4 with (low resolution, top) with run cbl6 (nigh 
resolution, bottom). 



